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Abstract 



In this paper we use microscopic arguments to derive a nonlinear 
Schrodinger equation for trapped Bose-condensed gases. This is made pos- 
sible by considering the equations of motion of various anomalous averages. 
The resulting equation explicitly includes the effect of repeated binary in- 
teractions (in particular ladders) between the atoms. Moreover, under the 
conditions that dressing of the intermediate states of a collision can be ig- 
nored, this equation is shown to reduce to the conventional Gross-Pitaevskii 
equation in the pseudopotential limit. Extending the treatment, we show first 
how the occupation of excited (bare particle) states affects the collisions, and 
thus obtain the many-body T-matrix approximation in a trap. In addition, 
we discuss how the bare particle many-body T-matrix gets dressed by mean 
fields due to condensed and excited atoms. We conclude that the most com- 
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monly used version of the Gross-Pitaevskii equation can only be put on a 
microscopic basis for a restrictive range of conditions. For partial condensa- 
tion, we need to take account of interactions between condensed and excited 
atoms, which, in a consistent formulation, should also be expressed in terms 
of the many-body T-matrix. This can be achieved by considering fluctuations 
around the condensate mean field beyond those included in the conventional 
finite temperature mean field, i.e. Hartree-Fock-Bogoliubov (HFB), theory. 



I. INTRODUCTION 

The observation of Bose-Einstein condensation (BEC) in alkalis 10-0 strongly moti- 
vates a description of the evolution of the condensate that takes full account of the micro- 
scopic nature of atomic interactions in a trap, both close to and far from equilibrium. The 
conventional description relies heavily on the well-known Gross-Pitaevskii Equation (GPE) 

also known as the nonlinear Schrodinger equation. In this equation, one assumes that 
the atoms are all effectively condensed and the atomic interactions can be accurately mod- 
eled by a pseudopotential, expressed in terms of the s-wave scattering length. This theory 
appears to make good predictions about the condensate's properties [f)Hl7H and it is obvi- 



ously desirable to put this phenomenological theory on a clear microscopic basis. In fact, in 
spite of its frequent use, a direct link of this effective interaction approach to microscopic 
properties of the gas, and a discussion of the effect of the mean fields on the intermediate 
states of a binary collision, appears to be lacking. In this paper we will address such issues 
in a derivation of a nonlinear Schrodinger equation based on these microscopic properties 
and what we believe are reasonable assumptions about them. In particular, we shall show 
how to take account of all possible repeated binary atomic collisions in the presence of a 
condensate, thus deriving a mean field evolution in terms of the many-body Transition (or 
simply T) matrix. 

Our approach is not limited to cases close to equilibrium, and therefore complements 
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earlier work in this area. For example, Bijlsma et al. |18[] have used a variational approach 
to calculate equilibrium properties, such as the normal and anomalous self-energies, in the 
many-body T-matrix approximation. In addition to this work, one of us has employed a 
functional formulation of the Keldysh theory to derive the time-dependent Landau-Ginzburg 
theory for the long-wavelength dynamics of an inhomogeneous weakly-interacting gas at 
nonzero temperatures, also in the many-body T-matrix approximation jl9 l. In this paper 



we will show, by different methods, how to derive a time-dependent generalized nonlinear 
Schrodinger equation]] in a trap for both zero and nonzero temperatures based on microscopic 
arguments. We shall explain that, to obtain this equation, we must neglect the effect of mean 
fields in the intermediate collisional states; we shall also discuss the physical significance of 
these effects. 

In an earlier publication p0| , we derived equations of motion for thermal averages of 



products of up to three (single-particle) fluctuation operators. In this way, we obtained a 
time-dependent version of the Hartree-Fock-Bogoliubov (HFB) equations p3|-p6|j in 



terms of actual interatomic potentials, and further generalized them by considering more 
complex anomalous averages (triplets). This closed system of equations can be used in 
two different ways: In the first approach, all averages of products of fluctuation operators 
evolve on similar timescales, so that the equations need to be solved self-consistently. Such a 
treatment enables us to investigate the possibility of further order parameters being present 
in our system. We point out that the possibility of pairing as a competing transition to BEC 
in the case of attractive interactions has already been investigated in |27j . The equations of 



20| further allow for the possibility of three atoms grouping together, such as for example, 
condensation of triplets p8| . However, there are also situations (in particular a dilute gas 
with repulsive interactions), in which the higher order (anomalous) correlations vary on 



1 To avoid confusion, we shall henceforth refer to the conventionally used equation as the GPE 
and the one we shall be deriving here as the NLSE. 
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faster timescales, enabling us to formally eliminate them from the equations of motion, i.e. 
by integrating over their effect during collisions. In this paper we will show how this formal 
(adiabatic) elimination of the pair correlation in our HFB equations gives rise to a nonlinear 
Schrodinger equation that includes the non-dressed repeated binary interactions (ladders) 
at zero temperature. We shall furthermore extend our treatment to nonzero temperatures, 
for which we shall also need to eliminate the triplets mentioned above. 

The equations of motion for averages of products of (up to three) single-particle fluctua- 
tion operators have been given in appendix A, although we refer the reader to [^] for more 
details. We emphasize that these equations do not merely bring the bare particle ladder 
interactions into our formalism, but also include more complex effects, such as dressing and 
damping of intermediate states during collisions. In this paper, we discuss how — and in 
what limits — these microscopic equations reduce to the conventionally used phenomeno- 
logical Gross-Pitaevskii expression. We also discuss the possibility of consistent theories 
outside these limits. Furthermore, we explicitly mention what processes must be neglected 
in order to obtain a NLSE. In fact, we shall see one has to neglect the effect of the condensate 
mean field and those due to the presence of excited atoms on the intermediate states of a 



binary collision P5fl . Identification of these terms shows how to explore the deviation from 
the bare-particle T-matrix due to dressing generated by the mean fields. This treatment 
enables us to make qualitative predictions about the validity regime of the phenomenological 
GPE, an issue we hope to address computationally for inhomogeneous gases in the future. 

One might expect the zero-temperature evolution of the mean field of a trapped Bose- 
Einstein condensate in the simplest version of the mean field theory, to be described — in 
the occupation number representation — by 

lh ^l = j2h^(t)z k (t). (i) 

Here z n (t) corresponds to the time-dependent mean value amplitude of the n th trap level 
which is obtained from the single-particle operators a n (t) according to the shift 



a n (t) = z n (t) + c n (t) . (2) 



Here the operators a n are defined by the usual decomposition of the Bose field operator 
$?(r,t) into any complete set of orthonormal single-particle states ip n (r), namely 

*(r,t)=$> n (r)a n (t) . (3) 

n 

In this simplest approach, the time-dependent condensate mean field depends on the Hartree- 
Fock hamiltonian h^°\t) given by 

h [ nl{t) = (n\E\k) +J2{ni\V\jk)z*(t)^(t) • (4) 

ij 

Here S contains the kinetic energy and trap potential and (ni\V\jk) represents the sym- 
metrized form of the actual (single-vertex) interatomic potential |3(J between a pair of 
colliding particles. This is defined in terms of ^ n (r) by 

(ni\V\jk) = l - {(ni\V\jk) + (ni\V\kj)} , (5) 

where 

(ni\V\jk) = J J drdrV;(r)^(r')y(r-r')V fe (r')^(r) (6) 

and V^(r— r') represents the actual interatomic potential experienced between two interacting 
atoms at each collisional vertex. Although Eq. (1) may on first sight appear to be equivalent 
to the conventional form of the GPE H] 

ih ^T-= (-^ + ^a P (r))$(r,t) + iVt/ |$(r ) t)| 2 <|.(r,t) , (7) 

this is not the case. The GPE takes account of all (repeated) collisional processes via a 
pseudopotential of the form V(r — r') = XJ^bix — r') where Uq = 4irh 2 a/m [31]]. On the 



contrary, the matrix element defined in (5) represents the instantaneous (i.e single-vertex) 
interaction between two atoms as depicted in Fig. 1(a). 

In this paper we will show that we can include the ladder diagrams of Fig. 1(b) into our 
formalism, by considering the evolution of correlations of products of fluctuation operators 
c n which are given in appendix A. These correlations influence the condensate mean field 



according to the exact equation |20| 



dz n 

ifc—^- = J2(n\E\k)z k + Y^(ni\V\jk) [z*ZjZ k + K jk z* + 2pjiZ k + Xijk] , (8) 

k ijk 

where we have defined the quantities 

Pji = (c\cj) , K jk = (c k Cj) = (CjC k ) , Xijk = (cUjCk) = (CkCjcj) . (9) 

In obtaining (8) we have merely assumed the existence of a mean field (in the sense of Eq. 
(2)) and that the dominant collisions in our atomic assembly occur pairwise. 

The main part of this paper deals with the microscopic derivation of a nonlinear 
Schrodinger equation for bare particles in a trap, and we shall discuss in what limits it 
reduces to simpler versions. The use of a NLSE necessitates a well-defined and thus slowly- 
varying condensate fraction, which in our formalism corresponds to a slowly-evolving mean 
field amplitude z n , in comparison with the duration of a typical binary interaction. This 
means that, in most of the analysis given below, we are limiting ourselves to the case of 
net repulsive interactions between the atoms, i.e. we are not allowing for BCS-type effects 
|27j to become important. When extending our treatment to nonzero temperatures, we will 



similarly need to assume slow evolution of the excited state populations, although this issue 
becomes rather subtle as we shall discuss in Sec. Ill C. In our subsequent calculations we 
shall discuss when the thermal correlations Kjk and Xijk vary on a much smaller timescale 
and show how they can be formally eliminated in our treatment. 

This paper is divided into five parts. In the first part of the paper we restrict ourselves to 
near zero temperatures, for which it appears reasonable to neglect all effects of excited states. 
We thus obtain a NLSE that includes the full (bare particle) ladder interactions, where 
this effective interaction potential arises from the adiabatic elimination of the anomalous 
correlation Kjk- This equation is then shown to reduce to the phenomenological Gross- 
Pitaevskii equation by means of the usual pseudopotential approximation. In Sec. Ill, we 
extend our treatment to nonzero temperatures. The first modification that comes about is 
the inclusion of the occupation of virtually excited states into the ladder diagrams. However, 
the finite number of excited atoms further necessitates the consideration of collisions between 
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condensed and excited atoms. We thus show how adiabatic elimination of 'upgrades' 
the condensate-excited state interaction potential to the ladder approximation. In Sec. 
IV, we show how the particle propagators get dressed by the condensate mean field and 
by interactions with other excited atoms and discuss the consequences of identifying these 
terms in our equations. We then discuss, in Sec. V, the validity of several mean field theories 
for both zero and nonzero temperatures and in particular point out the inconsistency of the 
conventional finite temperature mean field approach (i.e. what is commonly referred to as 
HFB). Finally, in Sec. VI, we discuss some further effects of the mean fields on collisions. 
The detailed form of the equations needed for the current discussion have been reproduced 
in appendix A, whereas appendix B gives an exact treatment of the condensate mean field 
evolution up to second order in the interaction potential. 

II. LADDER INTERACTIONS IN A BARE PARTICLE BASIS AT ZERO 

TEMPERATURES 

Let us initially assume that the atomic assembly is sufficiently dilute, so that we can 
treat the bare particles as weakly-interacting quantities. In this limit, bare particles form a 
reasonable basis set for the system and can thus be used for the description of the assembly. 
We shall take the non-interacting part S of the single bare-particle hamiltonian in a trap as 
diagonal, by writing 



If this basis accurately represents the state of the system (i.e. all dressing of states by the 
mean fields can be ignored), then the anomalous average n will evolve rapidly, thus enabling 
us to adiabatically eliminate it. The elements of p, however, vary on a slower timescale, as 
their evolution depends on an energy difference. This means that they should, in principle, 
be retained in the equation of motion (8) for the condensate evolution. Near T = 0, however, 
the occupation of excited states is extremely scarce, so that we can, in the first instance, 
neglect its effect on the condensate evolution. 




(10) 
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Within these approximations, Eq. (8) reduces to 

ih-Tf = huj n z n + ^2(ni\V\jk) [z*ZjZ k + n jk z*] . (11) 
dt ijk 

In order to adiabatically eliminate the pairing Kj k , we shall need its equation of motion 
given in appendix A. We note that in a bare particle basis at T = 0, Eq. (A6) essentially 
simplifies to 

where we have defined [2C] 

A kj = ^2(kj\V\ms) [z m z s + K ms ] . (13) 

ms 

We are now in a position to derive a NLSE valid at T = with the actual interatomic 
potential of Eq. (4) replaced by the two-body T-matrix. Such an equation will be shown 
to be equivalent to the conventional GPE if we replace the T-matrix by the zero-range 
pseudopotential. 

Integrating (12) and substituting into (11), we obtain 

dz - r dt' 

ijk 



th^ = hu; n z n + Y,(™\V\ms)z;z m z s + Y,(™\V\jk)z*(t) J -A e -K+^)(^') 

ims ijk ^ 

J2(jk\V\ms)\z m (t')z s (t')+K ms (t)\ , (14) 



x 

ms 



where the integration is to be carried out in the interval [to,t], in which the collision takes 
place. Since the gas is considered sufficiently dilute, the mean collisional duration is neg- 
ligible compared to the time between successive collisions, enabling us to take the limit 
to — ► — oo. In writing Eq. (14) we have made use of the approximation that the pairing k 
asymptotically relaxes to zero, so that Kjk{to) — » as to — oo. Physically, this is the same 
as saying that these anomalous interparticle correlations are dominantly created within a 
collision between condensate particles. 

A further crucial point to note is that, since we are trying to derive a NLSE, we are 
implicitly assuming that the mean value amplitude Zi(t) is a slowly varying quantity, that 



is, it remains practically constant over the time scale in which the correlation Kjk varies. 
Thus, we are justified in allowing the quantities Zj(t') present in the integrand of Eq. (14) 
to evolve freely with time by setting 

Zi (t') = z^e+^e-'') . (15) 

If we now intitially limit ourselves to second order in the interaction potential V, we note 
that we can disregard the K ms (t') term of (14), as its first contribution will be to order V 3 . 
Then, Eq. (14) reduces to 

dz 

ih-ij- = hu n z n + ^{ni\T 2B (E)\ms)z*z m z s , (16) 

ims 

which gives us the bare-particle two-body T-matrix in a trap to second order in the inter- 
action potential V, namely 

f 2B (E) = V + Y,V\jk) f —e-^+^-^-^ijklV . (17) 

jk 

Here E represents the energy of the colliding (bare) atoms prior to the collision in the 
centre-of-mass frame (i.e. in Eq. (16), E = % m + cu s )), whereas j and k correspond to 
intermediate trap eigenstates induced during the collision. However, we point out that it 
is precisely the presence of the term n ms (t') which generates T 2B to all orders in V, as can 
easily be seen by repeated use of equations (11)-(13). We have hence shown Eq. (16) to be 
valid to all orders in V, with T 2B satisfying the following Lippmann-Schwinger relation in a 
trap 

f 2B (E) = V + Y,V\jk) f —e^ + ^- E ^-'\jk\f 2B {E) . (18) 

jk 

This modified form of the NLSE accurately represents a collision between two condensate 
atoms occuring in vacuum, by taking account of all possible states via which the collision 
can proceed. Thus, (18) includes the ladder diagrams for bare particles as depicted dia- 
grammatically in Fig. 2(a). 

At this point, we could introduce the many-body T-matrix, by making use of the next set 
of terms in dn/dt (given in Eq. (A6)). However, before returning to this in Sec. Ill (A), let 
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us focus on the relation between this version of the microscopically derived zero-temperature 
NLSE (16) and the phenomenological GPE (7). 



A. Reduction to the Gross-Pitaevskii Equation 

Let us now show explicitly how we can obtain the GPE from our NLSE (16). To do this 
we shall need to approximate the (bare-particle) T-matrix defined in (18) in terms of actual 
interatomic potentials, by a phenomenological model, such as the pseudopotential approach 
discussed in detail in the classic paper by Huang and Yang |32fl . For the low condensation 
temperatures involved, we can limit ourselves to the lowest order contribution of the s-wave. 
We may thus write 

f 2B (r) ~ ^^5(r) , (19) 
m 

where a corresponds to the s-wave scattering length. 

Let us now use this approximation in (16) to explicitly reconstruct Eq. (7). In doing 
this, we shall also multiply the resulting equation by 4> n (r)/\fN and sum over the indices n 
(so as to obtain the desired $(r, t) on the LHS of our equation). We note that the factor of 
y/~N is needed so as to renormalize the wavef unctions, as the zi satisfy 

Y^z; Zt = N. (20) 

i 

After carrying out the integration with respect to the 5-function and rearranging, we thus 
obtain 



^JfM) = Jdv £<(r')0 n (r) S(r)J>(r 



Zk(t) 



+ NU J dr 



. N 



E<^(r)^E<Mr')^E<Mr)^§. (21) 



N j VN V VN 



Using the completeness relation for orthogonal states, this reduces to the GPE (7) for $(r, t). 

We shall now try to discuss the potential limitations of the GPE (7). To do this, let us 
briefly outline its conventional derivation, in which one starts with the Heisenberg equation 

11 



of motion for the Bose field operator *l/(r, t) in the limit of pairwise interactions. After using 
the appropriate commutation relations for the field operators and taking mean values, we 
are left with 

^^T 1 = (-^f + W')) *(r, f) (22) 
+ J drV(r-r')(tf + (r',*)#(r',t)$(r,*)) , 

where V(r — r) again corresponds to the actual potential experienced by a pair of interacting 
atoms (assuming a central potential). 

It is at this stage that the standard approximations are conventionally made (see e.g. 
|33[). One first uses the pseudopotential approximation V(r — r') = U S(r — v') and drops the 
interatomic correlations by setting (^(r, t)^(r', t)$(r, t)) = (^ t (r',t))(^(r',t))(^(r,t)). 
The two steps are interlinked and one could not carry them out independently as the resulting 
expressions would not be of any use. The second approximation is usually justified by 
arguing that the effect of all correlations thus neglected is taken into account in the effective 
interaction strength Uq. 

However, we stress that an effective interaction treatment can only be rigorously justified 
if the pseudopotential approximation is carried out after the two-body interaction potential 
has been 'upgraded' to the T-matrix, in the manner shown above. Although the resulting 
equation in this limit is identical to the GPE, we have given it a microscopic basis. Actually, 
we will use this basis to show that the terms ignored (in the GPE) are not always negligible. 

It is well known in the case of the homogeneous gas, that the use of a 5-function potential 
in a self-consistent treatment can lead to unphysical results, such as no depletion of the 



condensate [34]. This is due to the fact that such an approximation already contains an 



implicit assumption about the actual interaction operator, namely that it is acting on free 



particle states [[kJ . Put differently, the 5-function potential can only be used in combination 
with an ultra-violet cut-off in the theory. One must therefore be cautious of some predictions 
based on a zero-range effective interaction potential, unless this has been introduced in the 
rigorous manner on the two-body T-matrix (i.e. by means of (19)). 
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Furthermore, it is worth pointing out that the pseudopotential approximation (19) ex- 
plicitly ignores all momentum dependence, which appears in the problem as higher order 
terms in the s-wave scattering length a. This is an excellent approximation, valid, in the 
case of the homogeneous gas if a/A <C 1, where A represents the atomic de Broglie wave- 
length. However, we stress that our treatment is not limited to zero-range potentials and can 
handle spatially-dependent interaction effects that are not present in the phenomenological 
Gross-Pitaevskii model |36] . 



Thus far, in this section, we have only discussed the limitations of the GPE in terms of 
the treatment of the spatial and momentum dependence of the atomic interactions. Unfor- 
tunately, this is not the only problem we have to address. The main issue is that we have 
derived (16) based on a bare particle basis. This implies, for example, that we have ignored 
all dressing effects on the intermediate collisional states due to the condensate mean field. 
This obviously means that such effects are also implicitly ignored in the GPE, and implies 
that the GPE is strictly only valid if the atoms are colliding in vacuum. The reason why 
these mean field effects cannot be included in a phenomenological GPE based on the s-wave 
scattering length a, is that, in the experiments, a is spectroscopically determined in the 
absence of mean fields (i.e. effectively in vacuum). Our microscopic approach fully treats 
collisions in the presence of mean fields, and we believe it provides a natural formalism for 
obtaining expressions for both dressing and damping of the GPE 



III. LADDER DIAGRAMS AT FINITE TEMPERATURES IN A BARE 

PARTICLE BASIS 

Following the analysis of the previous section, it should be clear to the reader, that 
the GPE represents only a limiting case of the equations given in appendix A, since it 
neglects the effect of the medium in which the collisions occur (which becomes increasingly 
important as the atomic assembly becomes denser). Before we deal with these issues, we 
would, however, first like to extend Eq. (16) to finite temperatures, which corresponds to a 
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well-posed problem within our formalism. 

When extending our treatment to nonzero temperatures, a variety of new features arises: 
firstly, there will now be a non-negligible occupation of excited states, which leads to a 
modification in the scattering between two condensed atoms, as discussed in Sec. A below. 
At the same time, however, the evolution of the condensate mean field will also be affected 
by collisions between condensed and excited atoms, as will be dealt with in Sec. B. 

In this section we discuss the limit of weak interactions, in which the collisions can be 
described in terms of bare (i.e. non-dressed) atoms. The deviation from such a simplistic 
picture (due to dressing induced by mean fields) will be dealt with in Sec. IV. 



We shall now deal with the first of the finite temperature effects, namely how Bosc 
statistics affect the T-matrix, by taking account of the occupation of low-lying excitations 
during a collision between two condensate particles. This effect is expected to be so small 
at low temperatures, that we have chosen to neglect it altogether near T = 0. However, in 
order to take account of the occupation of virtual states accessed during a collision between 
two condensate particles, we need to modify Eq. (12) to 



We again formally integrate the above equation and adiabatically eliminate Kj^ from Eq. 
(11). Doing so we obtain, to second order in the interaction potential, the following equation 
of motion for the condensate 



A. The Many-Body T-Matrix for Weakly- interacting Particles 



(23) 



dt 



r 




+ J2(ni\V\jk)z*(t) f ^e-^ + ^ t -^Y.{3k\V\ms)z„Xt)z s {t) 




(24) 
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Here, the notation +{k <-> j} indicates the presence of an identical term upon interchanging 
the (intermediate propagation) labels j and k, and thus indicates the presence of both 
a direct and an exchange term (in much the same manner as those are included in the 
symmetrization of the interaction vertex (5)). Once again the k correlations ensure that 
this equation generalizes to all V. Hence, T 2B in Eq. (16) becomes replaced by the operator 
t, defined by 

t(E) = V + Y. V\i k ) j —e-^^-^-^ijkltiE) 
jk J 1 

+ En^>{ / —e-^ + ^- E ^-^p rr (t)(rk\ + (k ~j) \t(E) . (25) 

rjk I ^ J 

In the limit of weak interactions, we can treat the elements of p as diagonal, thus cor- 
responding to population of excited states. The off-diagonal elements will correspond to 
dressing of the bare atoms and will be discussed in Sec. IV. In this manner, we arrive at 
the Lippmann-Schwinger relation for the many-body T-matrix in a trap in terms of bare 
particles, in the form 

f MB (E) = V + J2 y\jk) J ^ e -*K+<^X'-'')[l + plit) + P l(t')) 

jk 

x(jk\f MB (E) . (26) 

This definition has also been depicted diagramatically in Figs. 2(b)-(c). The operator T MB 
of Eq. (26) is linked to T 2B of Eq. (18) by the equivalent definition 

f MB (E) = f 2B (E) + Y,T 2B (E)\jk) J ^ e -K+^-^)(*-*')[^.(t') + pl k (t)] 

jk 

x(jk\f MB (E) . (27) 

It is easy to see that, in the homogeneous limit, the expression (26) reduces to the 
well-known integral relation for the many-body T-matrix namely 

dk'' 



T MB (k', k, K; E) = V(k' -k) + J ^V{k' - k 



(2tt) £ 

l + N(f + k")+N(f-k", MB 



E - + zO 

m 
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T (k , k, K; E) . (28) 



Here k and k respectively correspond to ingoing and outcoming particle momenta, K to 
the centre-of-mass momentum and E represents the kinetic energy for the two atoms in the 
centre-of-mass frame. Furthermore, N represents the average occupation number of a single 
particle state (of specified momentum). 

B. Ladder Approximation to the Condensate-Excited State Interactions 

So far, we have only accounted for the occupation of states through which a condensate- 
condensate collision proceeds, in terms of the many-body T-matrix. The other important 
effect arising at finite temperatures, is the interaction between condensed and excited atoms. 
Obviously, such an effect must also be included in a NLSE describing the evolution of the 
condensate. To deal with this effect, we must now also consider the rest of Eq. (8) that 
has been ignored up to this point. To be more precise, we need to consider the effect of the 
additional contribution 

dz 

ih-£ = ... + ^2(ni\V\jk) [2 Pji z k + X ljk ] . (29) 

ijk 

Careful observation of the first few terms in the equation of motion for A^, Eq. (A7) 
suggests, by analogy to the arguments of Sec. II, that it is indeed the suitable quantity 
for upgrading the actual (single-vertex) interatomic potential for 2pjiZ k in Eq. (29) to the 
T-matrix level. However, we note that the triplet A is not included (i.e. A = 0) in the 
traditional mean field (i.e. HFB) theory. This shows that we need to go beyond HFB in 
order to rigorously obtain a NLSE for the evolution of the condensate mean field at non- 
zero temperatures, and we shall discuss this below. An extension beyond HFB is, strictly 
speaking, also essential for a consistent T = theory. 

We shall now also adiabatically eliminate []39] the quantity A^ appearing in Eq. (29). 



Consider initially the contribution 

ih-jifajk) = H^k + Uj - uJi)\ij k + ^(jk\V\ms) (2p mi z s + X irns ) . (30) 
Thus, Eq. (29) becomes, to second order in the interaction potential 
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<ft % = • • • + T,(™\V\ms)2 Pmi z s + £<m| V|jfc> / — e-^^-^-^ 

ims ijk ^ 

x Ys(jk\V\ms)2 Pmi {t)z s {t) . (31) 

ms 

As before, we assume that all mean value amplitudes Zi{t) vary on a much slower time-scale 
than the anomalous correlations. Furthermore, we limit our discussion here to the slowly- 
evolving diagonal elements p%5mi in a bare particle basis, for which p°(t') = p%{t). We 
thus obtain the desired upgrading of condensate-excited state interactions to the two-body 
T- matrix, valid to second order in V, just like in Eq. (17). Once again, we see that this 
expression can be generalized to all orders in the interaction potential, by taking account of 
the last term proportional to \ ims in the right-hand side of Eq. (30). 

Thus far, we have shown how to treat the condensate-condensate interactions to the 
bare-particle T MB level and the condensate-excited state contributions to T 2B , so that the 
condensate mean field evolves according to the equation 

th-^ = hu n z n + J2(m\T MB \ms)z*z m z s + 2j2(m\T 2B \ts)plz s . (32) 

ims is 

It is clear that a consistent finite temperature theory would require the latter contribution 
to be also expressed in terms of T MB . The subtle point that needs to be addressed here, 
is that such an expression will contain terms of order (p ) 2 ; hence, we expect (and we shall 
indeed confirm in our detailed treatment below) that in this situation, we should also take 
account of the scattering of excited states into the condensate, which modifies the evolution 
of the condensate mean value. 



C. Interactions Between Excited States 

In this section we shall discuss how to upgrade expression (32) to the many-body T- 
matrix for excited states. On first inspection, this appears to be quite straightforward, upon 
considering the first three terms of Eq. (A7), namely the contribution of 

ih^-(\ ijk ) = h((jj + uj k - uJi)\ ijk + ^2{jk\V\ms) {2p mi z s + \ ims ) 
at rms 
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+ \ ^{rk\V\ms)p jr [2p mi z s + X ims ] + (k <-> j) \ 



(33) 



to the equation for condensate evolution (29). 

However, we wish to point out that maintaining merely this contribution for A^-fc cannot 
be consistent. The reason for this is that the next term in Eq. (A7), namely 



- I ^2(ms\V\ir)p jm p ks z r . + (k j) \ (34) 

I. rms ) 



is an equally valid contribution to order V 2 in the expression for dz n /dt. Thus, a consis- 
tent approach must simultaneously take into account both contributions (33) and (34) in 
adiabatically eliminating A. The interpretation of the contribution (33) is straightforward: 
it leads to the renormalization of the condensate-excited state interaction to T MB . At first 
sight, it appears hard to interpret the physical significance of (34). To achieve this, we shall 
once again limit ourselves to the slowly-evolving excited state contributions p°. In this limit, 
maintaining both terms (33) and (34), we obtain for the contribution of excited states to 
the evolution of the mean value amplitude z n to second order in V, 

diZrr-) r * , . " . . f\ W * , . . . f (lit 



-i((j k +ujj~u)i-ujs)(t—t ) 



= ■ ■ ■ + 2j:(ni\V\is)plz s + 2j2{ni\V\jk) J -e 

is ijks 

Pi + P%P% + PiiPl ~ P'Ak] (t')(jk\V\is)z s (t) . (35) 



We can now easily identify the term in square brackets as 

{p% + 1)04 + i)p° - p%pUpI + 1) • (36) 

This contribution is known to give the correct amplitude for the scattering of quasiparticles, 
and can be used to predict condensate lifetimes. This expression is also in agreement with 



the results of the functional integral approach developed independently by one of us flT9| , f40" 

Hence, we can now identify the negative contribution in (34) as the term which ensures 
the correct factors for scattering into the condensate, due to the interaction of an excited 
state either with a condensed, or with another excited atom. Correct treatment of the 
factors of (36) to all orders in the interaction potential should lead to an additional term in 
the finite temperature NLSE, which should thus be expressable as 



= u n z n + Y.{nAT MB \ms)z:z m z s + 2 J>z|T MB |z S )p° , 



dt 

ims is 



ijks 

x(jk\f MB \is)z s . (37) 



Eq. (37) is depicted diagrammatically in Fig. 3. This equation is very important, as it 
includes the effect of kinetics due to interactions of excited atoms. At the moment it appears 
that our formalism does not generate the last many-body T-matrix, so that we merely obtain 
(ms\V\ik) . However, consideration of the derivation of the equations of appendix A reminds 
us that correlations of products of four fluctuation operators have only been treated by their 
mean field contributions, i.e. in terms of products of averages of two operators (or Wick's 



theorem, see Eq. (22) of reference p0|). This means that our formalism does not yet fully 
include effects on the population of excited states due to the interactions between two excited 
atoms. These can be taken into account by explicitly deriving the equation of motion for 
(c^c^cc) and carrying out the decomposition approximation on higher order correlations. We 



have indeed shown 0] that, in the limit of no condensation, such treatment gives rise to 
the well-known quantum Boltzmann equation JIT]. Thus, it appears reasonable that a full 
treatment of these more complex correlations might solve the apparent limitation of our 
formalism in terms of the derivation of Eq. (37). We shall return to this issue in detail 
in a following paper |37|| . However, since the last (kinetic) contribution arises due to the 
interaction of two excited atoms resulting in scattering into the condensate, this effect may 
actually be negligible in a sufficiently dilute and low temperature regime. 



IV. EFFECT OF MEAN FIELD ON COLLISION DYNAMICS 

All discussion so far has been in terms of weakly-interacting atoms, for which we can 
ignore the effect of the mean fields on the intermediate collisional states. In particular, 
in adiabatically eliminating the anomalous averages Kjk and Xijk, we have only considered 
certain of their contributions appearing in the equations of motion of appendix A. These 
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contributions correspond to those giving rise to T-matrices for both condensate-condensate 
and condensate-excited state interactions in the equation governing the evolution of the 
condensate mean-field. 

However, we know that the presence of mean fields (in terms of both condensate as 
well as excited states) will modify or 'dress' the intermediate collisional states. In this 
section, we shall limit ourselves to the regime where all effects of the triplets can be ignored 
(A = 7 = 0). This will be shown to correspond to the HFB description of the system 
in terms of actual interatomic potentials, with the effect of mean fields during a collision 
generating the quasiparticle (Bogoliubov) dressing. In Sec. V we will discuss the validity 
of the simple GPE, as well as that of other conventional dilute Bose gas theories and we 
shall argue that a consistent theory necessitates an extension of the conventional mean field 
theory (i.e. consideration of the triplets in our present language). 



A. The limit where A = 7 = 

Let us initially re-cast the equations of motion of appendix A in a simplified manner in 
the limit A = 7 = 0. In terms of the anomalous correlation Kjk, we note that Eq. (A6) 
contains the additional (dressing) contribution 



ih-^f = 1" H [VksKsj + Kksrfsj > (38) 

where rj is given by 

Vks = 2^2(kl\V\ts)(zlz t + p tl ) . (39) 
it 

From Eq. (A5), we also write down the equation of motion for p^. In Sec. Ill, we assumed 
k does not acquire a finite mean value in-between collisions (i.e. it relaxes to zero) and also 
ignored the off-diagonal p elements. To deal with the conventional mean field dressing, we 
shall now also allow for rapidly-varying off-diagonal p-elements Spji, i.e we shall substitute 
for pji the expression 
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Pji = l> [ i, (S j: + $Pji 



(40) 



We shall further assume that dp^/dt ~ 0. In these limits, we thus obtain the following 
set of equations for the rapid variations of 5p and 6k during a binary collisional process: 



iK- 



dt 



hu n z n + ^2(ni\V\jk) [z*ZjZ k + 2p° ii 8 i 

ijk 

+ ^2{ni\V\jk)z*5K jk 

ijk 

+2j2(ni\V\jk)z k 5 Pj 

ijk 



(41) 



ih—{5K kj ) 



fl(uj k + LUj)5K kj 

+ J2 [Vks$Ksj + SKksVsj 
s 

+ (l + P° kk + P° jj ) A kj 

+ [ S Pks ^sj + A ks 5p* sj 



(42) 
(43) 
(44) 



ihj t (5pji) 



h(ujj - LOi)8pji 

+ [VjrSpri - 5pj 
r 

+ (Pii - Pji) va 



(45) 
(46) 
(47) 



The physical significance of the various contributions appearing in these equations will be 
analyzed in detail below. First, however, we would like to point out that the limit A = 7 = 
of the equations of motion of appendix A corresponds to the time-dependent Hartree-Fock- 
Bogoliubov (HFB) equations, written in terms of actual interatomic potentials. These HFB 



equations can be alternatively generated by the following hamiltonian |26 
H Q = \ H{ h vi (cjA + c,cj) + (A M cJcJ + A; q c q c p ) } , 

where 



(48) 
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h pq = (p\h\q) = (j>\Z\q)+% g = (p|S|?) + 2'£<pl\V\tq)tfz t + p tl ) . (49) 

it 

This operator takes account of all possible quadratic terms in bare particle fluctuation 
operators q. In the literature, Hq is termed the quasiparticle hamiltonian ||26|| , as it describes 
fully an assembly of non-interacting quasiparticles. The dressing effect of the mean fields 
on the bare-particle states resulting from this hamiltonian is conventionally referred to as 
quasiparticle dressing, and we shall now discuss its consequences. 



1. Quasiparticle Dressing at Zero Temperature 

At T ~ there will be very few excited atoms in the assembly, so that we can approximate 
p\ ~ 0. We have already shown that the term ih(dnkj / dt) oc Akj of Eq. (43) corresponds to 
the two-body T-matrix being generated in the interaction between condensed atoms. The 
terms containing rj in both 5p and 8 k (i.e. contributions (42) and (45)) give rise to the 
dressing of Fig. 4(a). This represents the simplest dressing term due to the mean field of 
the condensate in the intermediate collisional steps, and corresponds algebraically to shifting 
and mixing of the intermediate states and frequencies. 

From Eqs. (44) and (47), we see that there further exist more complex contributions that 
depend on the coupling between bp and 8k. We expect these to correspond to the anoma- 
lous quasiparticle (Bogoliubov) dressing, in much the same way creation and annihilation 
operators get coupled via the quasiparticle transformation 



£ [uuk - v$] (50) 



where the b\ correspond to quasiparticle annihilation (creation) operators. 

Indeed, the combination of contributions (44) and (47) gives rise to intermediate dressing 
of the form of Fig. 4(c). We stress that the coupling between them implies that both 
diagrams in 4(c) are combined at any single intermediate step. To understand the form 
of the coupling, we have illustrated in Fig. 5 two typical diagrams in the evolution of the 
condensate mean field, along with their dressed equivalents due to the 8p-8n coupling. 
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2. Quasiparticle Dressing at Nonzero Temperatures 



Observation of Eqs. (42)-(43) and (45)-(46) shows the following new features arising at 
T > 0, where p% becoms non-negligible due to thermally excited atoms. The first one, which 
we have already discussed quite extensively, is the replacement of the two-body T-matrix by 
the many-body one, due to enhanced condensate-condensate scattering via occupied excited 
states. Furthermore, the 1] terms in (42) and (45) also allow modifications in the intermediate 
collisional states due to the mean field of excited states (in analogy to the condensate mean 
field of Fig. 4(a)), as shown in Fig. 4(b). 

The final modification we need to consider at finite temperatures is that due to the effect 
of the contribution (46). Adiabatic elimination of this contribution results in the generation 
of the so-called bubble diagrams shown in Fig. 4(d). These bubble diagrams can usually be 
neglected in dilute systems, but have recently been shown to be rather important near the 



critical temperature of the gas [43 



The interplay between all dressings of Figs. 4 (a)-(d) in the many-body T-matrix corre- 
sponds to what is conventionally refered to as the quasiparticle dressing for the interactions 
between condensed atoms. Indeed, the consistent combination of diagrams 4(a)-(c) with the 
T-matrix elements will generate the well-known normal and anomalous self-energies fo£n 
and hT,i 2 . However, we believe that the conventionally used finite temperature mean field 
theory does not actually include all these effects in a consistent fashion, as we shall argue 
below (Sec. V). 

Following the above discussion on the dressing of interactions due to the presence of 
mean fields, we shall now turn our attention to the analysis of conventional Bose gas theories 
appearing in the literature. 



V. CONVENTIONAL THEORIES AND THEIR VALIDITY 
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A. The Gross-Pitaevskii Equation 



In section IV, we discussed the effect of mean fields on the intermediate collisional states 
in the limit A = 7 = 0. This effect means we should use non-interacting quasiparticles, 
rather than bare particle states in our theoretical treatment. This dressing is not taken 
into account in the GPE, as was explicitly shown in its microscopic derivation in section 
II. However, this dressing can only be ignored in the limit of weak interactions, i.e. when 
nllo <C hu. In order to derive the GPE, we also have to ignore the effects of the occupation 
of low-lying excitations during a collision (many-body effects). These restrictions give us 
the sufficient condition for the validity of the GPE, namely 

nU , k B T < huj , (51) 

where hu is a typical energy separation between trap levels. 

Effectively, the GPE only treats condensate-condensate scattering in vacuum (thus ignor- 
ing both quasiparticle dressing and many-body effects). In this limit, the s-wave scattering 
length indeed contains all relevant collisional information. It is conventionally argued that 
the corrections to the GPE due to the presence of mean field are purely diluteness cor- 
rections. If the conditions of (51) are fulfilled the corrections are indeed small and this 
statement is close to the truth. However, even this statement is open to critisism [ft7|] . More 
generally, in the case of partial condensation (T 7^ 0), we maintain that the notion that 
diluteness is the issue is misleading. 



B. Zero- Temperature Bogoliubov-de Gennes Equations 

As mentioned above, the GPE takes no account of the finite T = depletion of the 
condensate, arising due to collisions between condensate atoms. This depletion is indeed 
very small in the case of typical experimentally-studied condensates, as has been shown 
explicitly, e.g. by Hutchinson et al. f45| . One can, therefore, calculate the frequencies of the 
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elementary excitations, to good accuracy, by finding normal modes of the linearized GPE, 
of the form 



$(r, t) = e"*"* [<fi(r) + u(r)e"^ + v(r)e iuJt ] . (52) 

Here /i corresponds to the chemical potential of the undisturbed ground state and the con- 
densate wavefunction has been represented by the condensate orbital </>(r). In addition, uj 
labels the frequency of the elementary excitations, whereas u(r) and v(r) correspond to the 
spatially dependent coefficients of the condensate's linear response to some driving field. 



This linear response approach has been discussed in JTI| , |15[ where it has been shown that 
the substitution (52) is completely equivalent to carrying out the Bogoliubov transformation 
|H on the fluctuating part of the Bose field operator ^ — (Jff), i.e. to the diagonalization of 



the binary-interaction hamiltonian for the assembly. 

Thus, one obtains a set of three static coupled equations for 0(r), u(r) and v(r), known 
as the Bogoliubov-de Gennes (BdG) equations. These equations have been used to predict 
condensate shapes, densities and the energies of elementary excitations at near zero tem- 



peratures. The results of the BdG equations appear to be in excellent agreement [TIJ with 
experiments as discussed in the general zero-temperature mean field theory review paper by 
Edwards et al. [15 |. 



C. The 'Essence' of Hartree-Fock-Bogoliubov 

We would like to point out that the equivalence between use of (52) and the diago- 



nalization discussed in ||11|| , [jl5|| is only true for a binary- interaction hamiltonian under 
the assumption V(r — r ) = UoS(r — r'). In this paper, we have argued that a ^-function 
pseudopotential approximation can only be imposed on an 'upgraded' effective interaction 
potential expressable in terms of the two-body T-matrix. This enormously complicates the 
discussion of the validity of certain approximate mean field theories, and generates confusion 
as to precisely what set of equations is implied by the term 'Hartree-Fock-Bogoliubov'. 
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We would argue, that the HFB equations are, in fact, defined in terms of actual inter- 
atomic potentials. Thus, we define HFB as the set of coupled z-p-K equations of appendix 
A (in which all interaction terms are left in terms of symmetrized matrix elements), in the 
limit A = 7 = 0. The time- dependent HFB equations defined in this way, can thus be 
derived from the quasiparticle hamiltonian (48) describing an assembly of non-interacting 
quasiparticles. We have already argued in Sec. Ill (A), however, that such a theory cannot 
consistently describe the interaction between an atom in the condensate with an excited 
one. The rigorous way of doing this, is by consideration of the A correlations. Bringing in 
the A correlations is equivalent to allowing the quasiparticles to interact with each other. 

We shall now discuss what people conventionally refer to as the finite temperature mean 
field, or HFB approximation. 



D. Finite Temperature Mean Field Theory and the Popov Approximation 



The mean field treatment extensively discussed in the literature avoids dealing with 
actual interatomic potentials. Instead, one conventionally expresses the above finite tem- 
perature HFB equations in terms of a 5-function approximation. In the static case, these 
equations then take the form given below. We would like to point out that these equations 
can be derived as a special case of our microscopic formalism if one includes triplets. The 
equations are 

r n 2 " 



2M 



V 2 + V t ^ p (r) -H + U A \^(r)\ z + 2n(r) ^(r) + U a m{v)^(v) = 0, (53) 



tuj{v) + Uq N {^(r)} 2 + m(r) Vj(r) = E jUj (r), 



(54) 



tvj{v) + Uq N {r(r)} 2 + rh*(r) Uj (r) = -E jVj (r), 



(55) 



where we have defined the following quantities 



= E{[M r )! 2 + MOI 2 ] NoiEj) + \v,(r)\ 2 } , 



(56) 
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m(r) = $>,(i>;(r) [2N (Ej) + 1] , (57) 

j 

and the operator 

t = "|^V 2 + K rap (r) -n + 2U Q [n |^(r)| 2 + n(r)] . (58) 
Here, No(Ej) = l/(er Ei — 1), and the above quantities are normalized in the usual manner 



TTH . These equations are being used at the moment to study the finite temperature excita- 



tions of Bose-condensed gases [Tiqjiql . Since these equations do not use actual interatomic 
potentials, we would rather not use the term HFB. In this paper, we shall call them the 
finite temperature Bogoliubov-de Gennes (BdG) equations, for want of a better name. The 
reason for making this distinction, arises from our view that the approximations made in 
obtaining Eqs. (53)- (55) actually contain more physics than the conventional mean field 
approximation. This can be seen from our microscopic approach: In Sec. Ill, we have ex- 
plicitly shown that the condensate-excited state interactions can be considered in terms of 
an effective T-matrix interaction only once the triplet A is taken into account. We stress 
that this T-matrix effective interaction is implicit in the above form of the finite temperature 
BdG equations (53)-(55), since they have been written in terms of zero-range potentials. It 
is clearly important to bear in mind that, when one includes the A terms, one also brings in 
other effects that have been ignored in the coupled BdG equations, such as further dressing 
of intermediate collisional states. The finite temperature BdG equations, therefore, have 
some inconsistencies built into them. How important these inconsistencies are in practice 
remains to be seen. 

Some discussion has been recently focused around the so-called 'Popov' approximation 
of these equations, which corresponds to setting m(r) = in Eqs. (53)-(55) p2| . This 
approximation has been used to compute the finite temperature excitations of a trapped 
Bose gas [[f5|,[|(J. In these papers, the equations being solved are referred to as the cou- 
pled HFB-Popov equations. We shall use the labelling Popov-BdG, since even the Popov 
approximation of these equations goes beyond the conventional mean field theory (due to 
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the implicit inclusion of the A correlations). In our treatment, we have obtained these equa- 
tions by adiabatically eliminating both anomalous averages k (or equivalently m(r)) and A 
(not explicitly present in Eqs. (53)- (55)). Thus, the Popov limit appears to be one way of 



obtaining a consistent theory. Indeed, Griffin [22] has shown that such a theory is gap less. 



However, using Eqs. (53)-(55) as they stand, with a nonzero value of m(r) is clearly incon- 
sistent, because we have shown that it is precisely the adiabatic elimination of m(r) which 
leads to the replacement of the interatomic potential by an effective interaction Uq. 

Let us now discuss the validity of the Popov approximation. The BdG equations have 
been written in terms of the quasiparticle coherence factors %(r) and Vj(r). This shows 
explicitly that they include the effect of the condensate mean field on the initial and final 
collisional states, which modifies the bare particle states into non-interacting Bogoliubov 
quasiparticles. However, the analysis of Sec. IV shows that in this limit, we must also 
consider other effects, like the dressing of intermediate collisional states and the occupation 
of these excitations. These many-body effects should be fully included in a consistent theory 
of condensate-condensate and condensate-excited state collisions. The Popov-BdG theory 
does not take account of these effects, and we would thus expect it to deviate from the 
actual description of the system. Indeed, a basis-set simulation of excitation frequencies 



using the Popov-BdG equations |46fl , has revealed large (and qualitative) differences from 
experimental data for temperatures above 0.6T C . We are currently working on determining 
a more appropriate set of equations — similar to what we have termed finite temperature 
BdG equations — that can be straightforwardly computed. 

Our treatment shows that many-body effects are not the only effects that are not included 
in the Popov-BdG description, as discussed below. 

VI. EFFECTS OF MEAN FIELDS BEYOND QUASIPARTICLE DRESSING 

The equations of appendix A, indicate that the dressing discussed in Sec. IV is not the 
only effect the mean fields have during a collision. The dressing of Sec. IV arises from 
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the quasiparticle hamiltonian (48), which replaces the interacting (bare) atoms by non- 
interacting quasiparticles. However, we believe the theory should allow the quasiparticles 
to interact weakly with each other. These interactions will lead to more complex dressing 
of the intermediate collisional states. However, this is not the only reason for going beyond 
a quasiparticle description. We have extensively shown in this paper that, to obtain the 
evolution of the condensate mean field at finite temperatures, we need to consider the triplet 
Xijk in Eq. (8). However, once including this, there appears to be no valid argument for 
neglecting the effects of the triplets in the remaining equations of appendix A. For example, 
Eq. (A6) shows that when adiabatically eliminating Kjk, we will also have to worry about 
triplet effects. 

We thus believe the triplets have two effects. The reader is by now familiar with the first 
one, which is the necessity of triplets for deriving finite temperature equations for condensate 
evolution. Secondly, to extend this argument even further, we also think that some of the 
triplets appearing in Eq. (A6) are actually needed in order to consistently combine the 
dressing diagrams of Fig. 4. We believe that a detailed study of the effects of the triplets 
may shed some light into why the dressed many-body T-matrix appears to go to zero as 
T — > JT8] , ^5| . This behavior implies a vanishing interaction (in the nonlinear Schrodinger 
equation) for collisions between two condensate atoms, which cannot be correct by itself, 
since we know that it is precisely these collisions which lead to the finite zero-temperature 
depletion of the condensate. We hope to discuss this, and related issues, in the future. 

VII. CONCLUSIONS 

In this paper we have carried out an in-depth analysis of our microscopic description of 
the behavior of Bose-condensed systems at finite temperatures. The equations of [|20[ have 
been re-cast in Appendix A, for the convenience of the reader. 

After a general introduction into our microscopic approach, we used the set of time- 
dependent Hartree-Fock-Bogoliubov (HFB) equations based on actual interatomic poten- 
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tials (i.e. without making any assumptions about them), to derive a nonlinear Schrodinger 
equation (16) for the description of the condensate mean field evolution at temperatures 
close to zero; in this limit, we can, to a good approximation, neglect the presence of excita- 
tions. We showed explicitly how to 'upgrade' the interatomic interaction potential (5) to the 
bare-particle two-body T-matrix (18) which includes all repeated binary collisional processes 
in vacuum. This is made possible by adiabatic elimination of the anomalous correlation (cc) 
which is not retained in the conventional derivation of the GPE. 

Furthermore, we explicitly showed in what limits the bare-particle NLSE (16) reduces to 
the GPE (7), thus providing a clear microscopic derivation. We then generalized this NLSE 
to finite temperatures by including the effects of excited states. This results in a variety 
of new features: Firstly, the occupation of excited states during a collision leads to the 
replacement of the two-body T-matrix for bare particles by the many-body one. Secondly, 
we must now also consider the effect of condensate-excited state interactions in the evolution 
of the condensate mean field. In order to do this consistently, we have shown the necessity 
to extend the conventional mean- field theory by explicitly including the triplet (c' cc) , which 
upgrades this interaction to the T-matrix. We have thus argued how to obtain a finite 
temperature equation for the evolution of the condensate mean-field, which includes kinetic 
contributions due to collisions between excited atoms. All such treatment is true in the limit 
of weak interactions, i.e. when nU <C fiuj. 

When this condition does not hold, we must take account of the effect the mean fields 
have on the intermediate collisional states. This leads to dressing of the states accessed 
during a binary collision. The presence of such dressing contributions makes it extremely 
difficult to adiabatically eliminate the anomalous correlations (cc) and (c^cc) in favour of 
a consistent equation valid in all limits. In fact, we point out that the triplet (c^cc) must 
be rigorously dealt with in a consistent T = theory, due to the finite (albeit negligible) 
condensate depletion. Our microscopic approach enables us to discuss the conventional 
mean-field theories currently used for describing the evolution of the condensate. In first 
instance, we have argued that the Gross-Pitaevskii equation is strictly only valid in the 
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regime nllo, ksT <C hu. Furthermore, we have argued that the conventionally used finite 
temperature Bogoliubov-de Gennes equations contain some inconsistencies if they are used 
beyond the Popov approximation. We have also argued that even the Popov limit of these 
equations (which is generally believed to form a consistent theory ||22j| ) may fail due to the 
neglect of many-body effects on the collisions in the gas. Indeed, recent simulations EH 
have shown huge discrepancies with the experiments at temperatures beyond 0.6T C . 

The analysis carried out in this paper shows precisely how hard it is to obtain a consistent 
mean field theory for the description of partially Bose-condensed systems. An alternative 
description that may not face the same difficulties is based on Popov's approach of describing 
the homogeneous Bose gas in terms of an effective condensate density and phase. Such an 
approach has already been discussed by Ilinski and Stepanenko |48| , and we hope more 
discussions in this area will appear in the future. 
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APPENDIX A: GENERALIZED MEAN FIELD EQUATIONS 



In this appendix we give the full set of self-consistent mean field equations derived in |]20 
subject to a decoupling approximation of correlations of four and five fluctuation operators. 
We have re-expressed these equations in a form that allows us to identify the physical impor- 
tance of all the individual contributions. We remind the reader of the following definitions 

Pji = , Kjk = (c k Cj) , X ijk = (cldjCk) , 7 ijfc = (CiCjCk) , (Al) 
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Vir = 2j2(jl\V\tr)(z:z t + p tl ) , (A2) 
it 

^kj = ^(kj\V\ms) [z rn z s + n ms ] . (A3) 

ms 

In writing down the equations of motion, we have assumed we are working in a bare particle 
basis, where the operator S = — (ft 2 V 2 )/(2m) + V tra p is diagonal, i.e (k\E\n) = hu n 5 nk . 
Furthermore, the equations given below are expressed in terms of the actual interatomic 
potential experienced between two atoms at each collisional vertex. We thus obtain the 
following set of equations 
dz n 

ih-^ = hu n z n + Y^(ni\V\jk)[z*z j z k + K jk z* + 2p ji z k + \ ijk ] , (A4) 

ijk 

ri PjrT]ri\ 



rms 

- E^l^l") ( 2 Ws* + Xjms Z r) , (A5) 



+A kj + ^ [p fcs A sj + A fcs p* 

s 

+ 2 {(&?"|y|ms) [2A rsj z m + imsjZ*] + (k ^j)} , (A6) 



+ ^2(jk\V\ms) (2p mi z s + \ ims ) 

ms 

+ E {(^|^|ms)p jr [2p mi z s + A ims ] + (fc <-> j)} 

rms 

- ^ {(ms|y|ir)p jm p fcs 2V + (fc <-> j)} 



32 



+ E { X ijsV* sk + \*iA* + (k <- j) } 

s 
s 

+ E {(kr\V\ms)2p mi (K js z* + X rjs ) + (k <-> j)} 

rms 

- 2(ms|Vlir) {^rk{2p jm z* + A* ms ) + (A; <-> j)} 

rms 

+ E {(^I^I^S) [2K mj (K* r Z s + \* sir ) + < r 7smj] + (k <-> j)} 
rms 

- ^(ms|V|ir) {2p jm A sfcr + (A; ^ j)} , (A7) 

rms 

at 

ms 

+ S {( ir l^l mS )/°J> [ 2 «fcm^« + 7fems] + (i <-> j)} 
rms 

+ ^(ir|y|ms) {/t mfc [/tjs< + 2A ris ] + (A; <-> j)} 

rms 

+ E 

s 

+(i ^j^k). (A8) 

(The term +(i <-> j <-> A;) in the last equation indicates summation of a// terms appearing 
in (A8), under cyclic rotation of the indices i, j, and k — except, of course, the first 'free 
evolution' contribution.) 

APPENDIX B: CONDENSATE EVOLUTION TO SECOND ORDER IN THE 

INTERACTION POTENTIAL 

In this appendix we give a systematic categorization of all contributions to the condensate 
mean field evolution, to second order in the interaction potential V. We shall carry this 
out by making use of the decomposition (40) of the single-particle correlation p into a 
slowly- varying diagonal (p°) and a rapidly-evolving off-diagonal element (5p). The implicit 
assumption made here, is that the atoms are weakly-interacting, so that the bare particle 
basis may be used, to a good approximation, for the description of the system. Furthermore, 
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we have already argued in Sec. Ill (C) that the equations of the appendix do not take 
account of the evolution of excited states due to quantum Boltzmann-type effects, which 
can be included by suitable treatment of the (cVcc) correlations. We shall therefore assume 
that during a collision dp%/dt ~ 0, which gives for the condensate mean field evolution to 
V 2 , the expression 

dz n 

ih—^ = hu n z n + Y]{ni\V\ms)z*z m z s (Bl) 
at 

ims 

+2j2(ni\V\is)p° ii z s (B2) 



IS 



1 ^ _ e -i(LO k +W j -LO m -UJ s )(t-t ) 



ijk ms J 1 

x{l + pl{t) +pl k (t))(jk\V\ms)z*(t)z m (t)z s (t) (B3) 

ijk s J 1 

x(l + P%{t) +pl(t)){jk\V\is) P Ut)z s (t) (B4) 

-EH^WE / — e - i (^+^-^-^)(*-*') 

ijfe s J 1 

x2p%(t)p kk (t)(jk\V\ l s)z s (t) (B5) 

+4^(m|T>|jA;)z fe (0E / — e - i ^+ w ™-^-^)(*-*') 

x - ^-(t')] (jm\ns) [z* m (t)z s (t) + , (B6) 

where t' acquires values within the range — oo to t. This equation has been diagrammatically 
depicted in Fig. 6, where we also give the contribution factors of each diagram. 

The terms in (B1)-(B2) represent the free evolution of the condensed particle, and the 
interaction of one condensed atom, either with another condensed atom (Fig. (a)), or with 
an excited one (Fig (b)). (B3) and (B4) give the lowest order corrections of the actual 
interatomic potential, due to its replacement by the many-body T-matrix, in both cases 
of condensate-condensate (Fig. 6(c)) and condensate-excited state interactions (Fig. 6(d)). 
The contribution (B5) — Fig. 6(e) — arises due to the interactions between two excited 
atoms and ensures the correct scattering amplitude factors for condensate evolution, as 
already explained in Sec. Ill (C). Finally, (B6) corresponds to the bubble diagrams of the 
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many-body formalism, which are shown in Fig. 6(f). 

At this point, we should really comment on the consistency of the above equation for 
condensate evolution. In obtaining this equation, we have not dealt with all quantities in 
the same manner. A fully consistent treatment would be to assign a slowly-varying and a 
rapidly evolving part to each of the quantities p, z and k (which would also generate two 
parts for both h and A), as well as the triplets A and 7 and then deal with them as a closed 
system. We shall return to this issue in a forthcoming publication |37 . 
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APPENDIX C: FIGURE CAPTIONS 



Figure 1: Fig. (a) depicts the contribution to the evolution of the condensate mean field 
(level n) due to the instantaneous interaction of two condensate atoms (labelled by j and 
k). In this, and all subsequent figures, the curly line represents the vertex between two 
interacting atoms located at r and r' respectively, and time runs vertically upwards. We 
have explicitly drawn the two contributions to the evolution of the condensate mean field 
due to the physical non-symmetrized matrix elements of Eq. (6) separately. Each of these 
diagrams is associated with a factor of (|) into our final equation for the evolution of z n , 
as defined by Eq. (5). In our notation, continuous lines with one 'free' end-point represent 
condensate particles, and have a factor z associated with them (the dashed condensate line 
gives no such contribution, as it is the one whose evolution we are monitoring in this paper). 
Fig. 1(b) shows the cumulative effect of repeated binary interactions of this type, or ladder 
diagrams, in which an arbitrary number p of loops (p > 0) may be present. Here, the final 
scattering into (condensate) states n and i has been mediated by other excited states (given 
by the vertical lines with arrows). 

Figure 2: Fig. (a) depicts diagrammatically the integral definition of the bare-particle 
two-body T-matrix, as defined by Eq. (18). In Figs, (b)-(c) we have illustrated the definition 
of the many-body T-matrix according to (26). The new feature introduced in (b) is that one 
now takes account of the occupation of the intermediate states through which the collision 
proceeds. Thus, the thick arrows indicate that the collisions actually occur in an atomic 
medium, as opposed to the vacuum. In particular, Fig. (c) indicates explicitly how this 
excited state occupation is taken into account according to the many-body factor (1 + p° ■ + 
p2 fe ), with the p° term indicated by the double arrow. 

Figure 3: This figure represents diagrammatically the evolution of the condensate mean 
value amplitude z n as given by Eq. (37). We have explicitly illustrated diagrams correspond- 
ing to direct and exchange terms, as well as the associated contribution pre-factors. These 
factors have been suppressed in our mathematical analysis, by means of the definition (5) 



40 



of a symmetrized matrix element. The factor p° indicates the population of excited states. 

Figure 4: This figure illustrates all different quasiparticle dressing effects on the unoccu- 
pied intermediate propagators of the two-body T-matrix (Fig. 2(a)). Such dressing during 
an atomic collision arises due to the effect of (a) the condensate mean field, or (b) the mean 
field of thermally excited states. Additionally the 5p-5n coupling generates the anomalous 
terms of Fig. (c), corresponding to the creation, or annihilation of two condensed atoms, in 
favour of excited ones; these contributions do not appear separately, but are complementary 
of each other. Fig. (d) indicates the extra dressing due to the bubble diagrams (46) which 
are associated with a factor (p° — We note that all above dressing effects could also be 
very straigthforwardly included in the many-body loops of Fig. 2(b)-(c). 

Figure 5: This figure illustrates the cooperative action of the 8p-5n coupling represented 
in Fig. 4(c), by means of two typical terms in the evolution of the condensate mean field, 
and their corresponding dressing due to this coupling: (a) shows the effect on the scattering 
between a condensed and an excited atom, whereas (b) corresponds to the dressing on a 
typical two-body ladder diagram. 

Figure 6: This figure illustrates the evolution of the condensate mean field to second 
order in the interaction potential. Figs, (a)-(d) show the many-body T-matrix terms in the 
interactions (to order V 2 ) between two condensed atoms (Figs. (a),(c)), or one condensed 
atom with an excited one (Figs. (b),(d)). Fig. 6(e) shows the contribution to condensate 
evolution, due to the interaction of two excited atoms. Fig. 6(f) corresponds to the bubble 
diagrams, which have a factor — p®^ associated with them (indicated by the hollow 
arrows). Continuous lines with one 'free' end-point have a factor z associated with them, 
and the thick arrows due to excited state occupation have been defined in Fig. 2(c). 
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